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- - - Abstract 
\Q ■ 

GN ! In this paper we solve the following problems: (i) find two differential operators P and Q satisfying 

[P,Q] = P, where P flows according to the KP hierarchy dP/dt n = [(P"/ p )+,P], withp := ordP > 2; 
(ii) find a matrix integral representation for the associated r-function. First we construct an infinite 
dimensional space W = spa,n c {'ipo(z),ipi(z), . . of functions of z £ C invariant under the action 
of two operators, multiplication by z p and A c := zd/dz — z + c. This requirement is satisfied, for 
arbitrary p, if ipo is a certain function generalizing the classical Hankel function (for p = 2); our 
representation of the generalized Hankel function as a double Laplace transform of a simple function, 
which was unknown even for the p = 2 case, enables us to represent the r-function associated with 
the KP time evolution of the space W as a "double matrix Laplace transform" in two different ways. 
One representation involves an integration over the space of matrices whose spectrum belongs to a 
wedge-shaped contour 7 := 7+ + 7~ C C defined by ^ = M + e ±7 "/ p . The new integrals above relate 
to the matrix Laplace transforms, in contrast with the matrix Fourier transforms, which generalize 
the Kontsevich integrals and solve the operator equation [P, Q] = 1. 
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It is a long-standing puzzle in the theory of 2(i-gravity to find an adequate description of gravitational 
coupling of (p, q) minimal models. One part of it is to find two differential operators P and Q of 
order p and q respectively, such that [P, Q] = f(P) for some function /. In the simplest case of q = 1 
' and / = 1, such description is provided by 1-matrix models, especially by the Kontsevich integral 

and their generalizations; see |l], |f9 , [25 1 . Going along the chain, 2d-gravity — > equilateral triangles 



discrete matrix models — > Kontsevich models, this approach has lead to the discovery of integrable 
structures for non-perturbative partition functions, which take the form of r-functions of the KP 
hierarchy (see J?], |25|, |3f[] for review and references). While similar results are believed to be true in 
the general (p, g)-case, the Kontsevich integral counterparts are still unknown. Note that a minor 
modification of the generalized Kontsevich integral can be interpreted as a duality transformation 
between (p, q) and (g,p)-models Q. 

So far the most promising approach for finding integrable structures in the general (p, g)-case seems 
to be the one initiated by Kac-Schwarz in the case q = I and f — 1. So, the general problem comes 
in two stages: (I) Find a point in Sato's Grassmannian invariant under two symmetry operators, 
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satisfying some commutation relation; the existence of such a plane leads to a system of differential 
equations specifying the wave function "J and thus to an algebra of constraints for the r-function. 
(2) Find a matrix integral representation for this r-function. Note a matrix representation, beyond 
the case q = 1 and / = 1, if it exists at all, was unknown. 

The purpose of this paper is to find a r-function and a matrix integral representation for the 
equation [P, Q] = P for q = 1 and arbitrary p . Remarkably, the matrix integral representation can 
still be found, but it is far less straightforward and considerably more involved, than the ordinary 
Kontsevich integral. 

The message is the following: whereas the case [P, Q] = 1 is described by general matrix Fourier 
transforms, a solution to [P, Q] = P is related to double Laplace transforms. While it is not known 
whether this solution has immediate physical relevance, it may help to shed some light on the (p, q)- 
case and on the matrix representations of the corresponding r-functions. In particular, what are the 
proper multimatrix generalizations of the Kontsevich integrals? 

Note this problem has come up in the physical literature, in various different contexts: unitary 
matrix models have been written down, leading to equations [P, Q] = P for differential operators 
P and Q in the double scaling limit; see the studies of Dalley, Johnson, Periwal, Minahan, Morris, 
Shevitz, and Watterstam [Q, ||, |2?], |28|, |2^, |23||). In the mathematical context (inverse scattering and 
monodromy preserving transformations), see Ablowitz, Flaschka, Fokas and Newell JD], ||, [hJ). The 
solution provided in our paper is new and does not require any scaling limit. 

Consider the problem of finding a differential operator P of order p and another differential 
operator Q satisfying 

[P,Q} = f(P), with ^= f(z) e C[z] . (1) 

When P is (formally) deformed with respect to the KP flows, i.e., dP/dt n = [(P"/ p ) + , P], one can 
introduce the corresponding deformation of Q which preserves Eq. (|l|). Hence ([!]) can be considered 
as a condition on a solution of the p-reduced KP hierarchy. 
The basic ingredients of this construction areQ 

o j4:C((z -1 )) — » C((z _1 )) which increases the order of an element of C((z~ 1 )) in z exactly by 
one, 

so that W := span c {^o> Aipo, ^4 2, 0o, ■ ■ ■} belongs to the big stratum of the Sato Grassmannian and 
satisfies AW C W, such that 

o -0Q satisfies the differential equation v(z)tpo = F(A)ipQ for some v(z) £ C((z -1 )) and F(Z) G 
C[Z], so that v(z)W C W also holds. 

Let ^ be the KP wave function corresponding to W. The above conditions lead to the existence 
of differential operators Q and P in x such that Q"i> = and P^P = v(z) x S. If A coincides with 
d/dv = (l/v')d/dz up to the conjugation by a function, then we have [P,Q] = 1. And if ipo is 
defined by a Fourier transform and the action of A on it can be expressed in a suitable way, then 
the corresponding Hermitian matrix Fourier transform, properly normalized, is the corresponding 
r-function. See Sect. || for details. 

The matrix integral approach to (|l|) has so far needed ord Q — 1 at the initial point of the formal 
KP time flows, requiring deg z f(z) < 1. The degree case can be reduced to [P,Q] = 1. In this 
paper, we provide a solution to the degree 1 case, or the next simplest instance of (|l[), which can 
clearly be reduced to 

[P,Q]=P, (2) 

with differential operators P and Q. As in the case of [P, Q] = 1, we write the r-function of its formal 
KP deformation explicitly in terms of a matrix integral. 

1 C[[x]] := {Yln°=o a " x " I «n S C} is the ring of formal power series in x, and C((x)) := {J] _ 00<Kn<00 a„x n | a„ <= 
C} is the ring of formal Laurent series in x. 
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Definition 1 Let — 1 < c < 0, p £ Z, p > 2. Let W &e i/ie linear span 

W = spa,n c {i/jo(z),tpi(z),tp 2 (z), . . .} , 
o/ generalized Hankel functions, 

pC+l poo z -cf uz \k e -{v.-l)z 

^ k ^ = T{ZcjJ 1 („» _ i)c+i du ' fc = 0,1,2,..., (3) 

a/so representable as double Laplace transforms 

c+l roc poo 

Mz) = ^n ziP ~ 1){C+1)e J dxx ° e ~ XZ J dyf k (y)e~ x y p (4) 

of the functions 

fk(y) = (<: k+1 e- Cy -C k - 1 e- rly )v k , k = 0,1,2,..., urift C:=e**/p. (5) 

Using the asymptotic expansion tpk(z) = z k (l + 0(l/z)) 6 C((z -1 )) as 5ftz — > oo, W defines a point 
of the Sato Grassmannian Gr. Let $ and r be the wave (formal Baker-Akhiezer) function and r- 
function, respectively, associated with the KP time evolution VV* = tiZl >V; see Sects. [j] and ^. 
Then we have 

Theorem 1 

*{x,0,z) = e B *il>o({l-x)z), (6) 

and it satisfies 

L\x-1, - z p ^j V(x, 0, z) = and U, - lj - (x - l) p ^j V(x, 0, z) = , (7) 
where L(z, d/dz) is the monic differential operator 

K'-s) »?(n(^+H-'' , 2('s + '- < )) - (I" 

Note that for p = 2, d/dz) = (d/dz) 2 - (c 2 + c)/z 2 . 

Theorem 2 Let 7ijv &e tfte space of N x N Hermitian matrices, and TC^ the subspace of Hn of 
positive definite Hermitian matrices. The corresponding r-function evaluated at 

t n := tr Z~ n , for n — 1,2, ... , and with an N x N diagonal Z, 

n 

is given by the following (normalized) double matrix Laplace transform: 

L+ dXdetX c e- trZPX L+ dYS (y)e-^ XYP 

T(t) = Sl(t)^£ 7^ r , 

J HN dX exptr p+ { j a 

where ( )2 denotes the terms quadratic in X , 

So(y) ■= ||^det(/ fc _ 1 ( 2/l )) 1 < 2!fc < JV and := det(Z^ c+1 ^)e tr z , 

where y = (y x , ...,y N ) are the eigenvalues of Y, y p = (y[, y p N ), and A(y) := Y[i>j{Vi ~ Vi) = 
det(y^ 1 ) ij , and f k -\ are as in 
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The function r(t) also has the following matrix integral representation 

J Hl m{dW) f n+ dX del X c (A{ W p)/A(w)) e ^(z~w) e tr x(W-z") 

T ~ f ~iv Tl (X + Z)p+ 1 \ ' 

J^rfXexptr(- ^ t_2 ^ 

integrated over the space of matrices 

H 7 N = {W= UDyU' 1 | 17 G XS{N),D 1 := diag(tui, . . . , w N ) G (7)^}, 

where 7 denotes a wedge-shaped contour in C, defined in Sect. Q ('see i* 1 117. [i[) ; m terms of a complex- 
valued measure 

m{dW) = dUdw J] (w, - Wj f. 

l<i<j<N 

Theorem 3 (i) The algebra of stabilizers ofW, 

S w ■= {<f>( z ,d/dz) G C({z~ l ))[8/dz] such that cf>W C W}, 
is generated by A c := z-S^ — z + c, z p and £ := z~ p F(A c ), where F(u) = 1 (w — i) — cpYlo~ 2 ( u ~~ *)■' 

5 w =C[yl c ,zP,e]cC((z" 1 ))[a/9z]. 
Moreover, W = C[Ac]^0) <™<^ V'o satisfies the differential equation 

F(^o = H%(z). (8) 

fn) .A family of solutions to the operator equation [P, Q] — P is given by the differential operators P 
and Q, defined equivalently by 

Pt> = z p y, Q® = -A C V, (9) 
P 

or by 

P = s(4-) S- 1 and Q = -(MpV p -P 1 / p + c), 
\dx J p 

where M = S kik(d/dx) k ~ 1 ) S~ 1 , tk — tk + &k,i%, with wave operator 



S = 



r(i- [{d/dx)-i}) 



r(t) 

(Hi) The function r(t) satisfies, in terms of the W -generators in Eq. $2(\), the following constraints 
]T a m ,i P.) ^ " W^_. T (t) = a m ,», e r(t) , to, n = 0, 1, 2, . . . , (10) 

0<i<m ^ 
0<j<« 

/or some constants a mt n iC , where the constants a ni i are defined by the formula^ (x ■ d/dx) n = 
'Yli i=0 Oi n ,iX l (d/ dx) 1 . In particular, setting m = 1, r(t) satisfies Virasoro constraints of the form 
(withW$=Z i+j=nr) :jPjP:) 

1 9 \ 

o W »p-«i «i,n,c r = 0, n = 0,l,2,.... (11) 

Remark 1 The constants a mi „ iC m can all be calculated; in particular, the Virasoro constraint 



(11) for n = becomes: 



» a a C (i + c )(p-i) \ 
V at* ^1 2 ) r ~ u - 



More explicitly, a n ^ = Jj- X]}=o (*) ( — -0* Note that it vanishes if n > Q and i = 0. 
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5.2 Proof of Theorem 2 



5 Proof of the Main Statements 

5.1 Proof of Theorems 3 and |l] and Remark [j] 



1 The KP Hierarchy 

Throughout, i is a formal scalar variable near 0, and z is a formal scalar variable near oo. If 
g(z) = cz q (l + 0(z -1 )), c 7^ 0, then ovd z g{z) := q is the order of g(z). 

Throughout, we denote d/dx by D. The algebra of ordinary pseudodifferential operators in x is 
denoted by T> (the word "in x" may be dropped if there is no fear of confusion), with its splitting 
V = V + +T>- into the subalgebras of ordinary differential operators and of ordinary pseudodifferential 
operators of negative order: 



-oo<i<n 



n € Z arbitrary, a, = Oi(x) 



A = a i Di eV A + = X! ajjD1 e V + and ^" = ^ ~ A + e T> - ■ 

The ring T> acts on the space of functions of the form X)-oo<i<oo o J i{x)z l e xz simply by extending the 
formulas D n e xz = z n e xz and A{Be xz ) = {Ao B)e xz , A, B e V. When A e V+, this definition of 
A(Be xz ) coincides with the usual action of A, as a differential operator, on Be xz as a formal series 
in x with z-dependent coefficients. 

A pseudodifferential operator in x may depend on the KP time variables t — (t\, t%, . . .) introduced 
below, but not on z unless otherwise noted. We are not specific about the regularity of the coefficients 
of pseudodifferential operators. The operators S, L, M etc., associated to a point W of the big 
stratum Gr° of the Sato Grassmannian (see below) have regular (i.e., formal power series) coefficients; 
otherwise, the singularities of those operators can be controled by the Schubert stratum to which 
W € Gr belongs. In particular, there exist n, m > such that x n S and S~ 1 x m at t = have regular 
coefficients. See 0] for details. 

As in 0), we set t = (x + tx,t2, £3, • ■ ■), and 

- ( d_ ld_ ld_ 

~ \dh' 2 dt 2 ' 3 dt 3 ' ' 
The elementary Schur functions p n are defined by exp = So° Pn(t)z n . 

1.1 KP hierarchy 

The operator L = L(t) = D + J2J=-oo o,j(x,t)D^ e T>, with t = (ti,t 2 , ■ ■ .), subjected to the KP 
equations 

?- = [{L n ) + ,L], n = l,2,..., 
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is known to have the following representation in terms of an operator S G 1 + T>- called the wave 
operator, and the associated, formally infinite order pseudodifferential operator 

W := Se^^ UD \ 

as follows: 

L = SDS- 1 = WDW-\ (12) 
OS dW 

m;r- {Ln) - s > and ou: = {Ln)+w - 

The wave function 

:= V(x,t,z) := We xz = StP^*"' , (13) 

where U = + Si^x, satisfies 

d<Z> 

L* = z* and — = (£")+*, (14) 

and has the following representation in terms of a scalar-valued function associated to S called the 
tau function r: 

*(«,*) - r^Lplesrw 

r(t) 
^ rftl 



n=0 



implying in view of (II 



r(t - [-D -1 ]) ._ g Pn(-d)r(i) D _ n _ (15) 



Moreover, using (|l3|), we have 



= ^-We** = W-H-e" = Wxe xz = WxW' 1 ^ , 
oz oz Oz 

thus leading to the operator 

M := WzW- 1 =SeV t **xe-V**° h S- 1 = s(x + ^kb k D k -^S- 1 

= S (j2kt k D k -^J S- 1 (16) 

satisfying 

M# = (d/dz)$ and [L, M] = W[D, x^' 1 = 1 , 
and for any formal series / = 

f(M,L) = Wf(x,D)W- 1 . (17) 
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1.2 Symmetries 

Consider the Lie algebra Woo of operators 



Woo ■= C[z, z 1 ][d/dz]=span c iz a \JI^J 



a, [3 S Z,0 > 



and its completion Woo := C((z 1 ))[d/dz] in the z ^-adic topology, for the customary commutation 
relation [ , ]. Acting on 'J, we have 

z a (d/dzfy = M p L a q> , (18) 
motivating the definition of the following vector fields, called symmetries, on "J: 

Y za(d/dz) ^:= (M?L a )^. 

We require that these flows act trivially on parameters x, t, and hence on S' _1 MS' = kikD k ^ 1 , for 
instance. 

Lemma 1 There is an injective homomorphism of Lie algebras 

{Lie algebra of vector fields 
on the manifold of wave functions ^ 
commuting with the KP flows d/dt n 

* a (^) — Y,„ (9/az)3 * = (M^)_*, 

^z"(d/dz)ff, ^z^'(d/dz)f = ^[z^(d/dz)^,z"'(d/dz)0'] • 

This definition differs from the one in Q by the sign. Here this definition is chosen to make it 
consistent with the natural action of mJoo on the Grassmannian discussed in the next section, rather 
than its negative. These vector fields induce vector fields on S and L = SDS^ 1 , as 

Y za(d/dz)0 (S) = (M L a ).S 

and 

Y z „ (a/az)3 (L) = [(M^ a )_,L]. 

Proposition 1 ([§]) We have 

, n,<eZ,B>0, (19) 



(M"L"+')-VE = _ 1} ^TWf^W 



ii — >ti+x 



where the Wi" , £/ie generators of the Woo-algebra, are the coefficients in the expansion of the 
vertex operator 

( 00 \ ( 00 \~ « - t o \ 
X(*,A,aO := exp - A*)*, exp £ — ( 1 



91, 



E^fc! E A- £ - fc < } , mtt wj 0) =^ )0 . (20) 

fc=0 ' fc-oo 
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2 Grassmannian 



Let H := C((z" 1 )), H+ := C[z], and H- := z^C^' 1 ]], so that H = H + © H- . We denote by Gr 
the Grassmannian manifold of linear subspaces W of H of relative dimension with respect to H + , 
i.e., the natural map 

7r w : W <=» H -2U ~ tf + 

being Fredholm of index 0. Gr° := {W £ Gr | 7ryy is isomorphism} is the big (open) Schubert 
stratum of Gr. 

Given a wave function = ^(x, t, z), let W be the point of Gr defined byQ 
W - sp an c |^— ^(0,0,z) 



^JiH kjN 



J = 0,1, 2,. 

N>0, ./; ./ V -0 



The first line guarantees W S Gr, and the second line follows from the first by using the second 
equation in (^), i.e., the KP time evolutions of 'J'. Hence up to the i-adic completion we have 

dY , 

.7=0,1,2, 



W = spanc|l — 1 *(0,M 

so that, letting ip = e~^< tiZ *v]/ and 

w t ._ e - £ uz' w = sp an c { (d/dx) j V>(0, i, z) | j = 0, 1, 2, . . .} , 

we have ?/> = (7Twt) _1 (l), i.e., "0 is the preimage of 1 by the map iry^t: W* — ► H + . 
The corresponding r-function r(t) is the determinant of the composite map 

w w* ff/fr_ ~ h+ , (21) 

where g denotes the multiplication by e - ^ tiZ . Given W, the determinant is well-defined up to a 
constant which is determined by the choice of a basis {0/c}^L O ' V'fe = zfe (l + C(z -1 )) for fc ^> 0, of 
W. We take {z fc }^o as ^ ne basis of More specifically, r(t) is defined as the limit as n — > oo of 
the determinant of 

W„ W -» #+ -» H+/z n H+ , (22) 

where the middle arrow is the composite map in (21), W„ = span c {0fe}j!~p , and the determinant 
is computed with respect to the bases {ipk}kZo °f W n and {z fc }]!~Q of H + /z n H + . The limit exists 
in the i-adic topology of C[[£]], i.e., for any multi-index a, there exists a positive integer n a such 
that, if n > n a , then the coefficient of t a in the determinant of (22) is independent of n, and gives 
the coefficient of t a in r(t). This finiteness property is an immediate consequence of the fact that, 
expanding r(t) in terms of Schur functions, the coefficients give the Pliicker coordinates of VV. See 
H for details. 

The wJoo-action on \& becomes the natural action of Woo on Gr: As an ordinary differential operator 
in z, each A £ Woo acts on H, which defines a vector field on Gr. 



2.1 Stabilizers 

Given W £ Gr, we shall call 

S W := {Q ■= Q(z, d/dz) £ Woo | QW C W} 

the stabilizer of VV. In this subsection we shall observe basic properties of the stabilizer which can 
be obtained without referring to matrix integrals. 

3 If \I/ is singular at (x,t) = 0, we need to replace (9/9x) J \E'(0, 0, z) in the first line by (d/dx) J I x n ^(x, 0, z) ) 

r— I V / \x=0 

for some n > 0, and make a similar replacement in the second line (see [M for details). We chose to write the formulas 
for VV G Gr° for simplicity. 
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Lemma 2 LetW G Gr and A := E-oo<i<oo,o<j«oo CijZ^d/dzY G Woo- 2/ 

iWcW, (23) 

f/ien 



g A := 51 cyM'L* G 2?+ 



0<j<cx) 

Conversely, if Q G 2?+ is of i/ws form, i.e., Q — Qa for some A G w M , i/ien i/izs A satisfies ((£ 
Proof: We have 

A*(t, «) = (24) 

by definition. Since AW C W, and since the Taylor coefHcients (or Laurent coefficients if W £ Gr°) 
in x of * generates W, A* is a C[[x, t]] -linear combination of D^>, Z? 2 *, . . . , i.e., A* = for 
some Q G 2?+. Hence, since (^) determines Qa uniquely, Qa itself must be in 2?+. Conversely, 
suppose Qa G T>+, and let 0, z) = E f n (z)x n be the Taylor (or Laurent) expansion of \&(ar, 0, z) 
at x = 0. Then each Taylor coefficient in x of Qa^ is a linear combination of {f n (z)}, and hence it 
belongs to W, so that by (^3) Af n G W for every n (the action of A on /„ is well-defined since A is 
a differential operator in z). Since {/„} is a basis of W, we have AW C W. ■ 

Corollary 1 Let p 7^ be an integer, and let Q G T> + such that a.d(L p ) N Q — for N ^> 0. TTien 
Q = Qa for some A G Woe swc/i iftai AW C W holds. In particular, a solution to the string equation 
always comes from a pair of A G uJqc and W G Gr, such that AW C W ( and z p W C W due to 
the extra assumption P = LP G 2?+ ^. 

Proo/: Writing Q = £ y Cij M^L\ let A = £ y c^z^d/dzY '. Since ad(L p ) Ar Q = we have ad(z p ) Ar A = 
0, which implies that A is a differential operator in z. Hence the "converse" part of Lemma |^ ap- 
plies. I 

Lemma 3 Let A, B G Wqo, tpo = 1 + 0(z~ x ) G 1 + 2T- and W G Gr. Suppose A acts on the 
monomials z k , k G Z, as 

Af fc = z fc + 1 (c fc + 0(z- 1 )) 1 
and Cfc ^ «/ fc > 0. T/ien i/ie following conditions are equivalent: 

(i) tpo G W, AW C W and SW C W; 

('ijj W = spancl^o, Ai/^o, A 2 ^0j ■ ■ ■}> fl^rf "00 satisfies the differential equations 

BA n ipo = F n {A)ipQ , n = 0,l,... (25) 

/or some F n (s) G C[s]. 

in particular, under these conditions W belongs to the big stratum Gr° of Gr. //, moreover, A and 
B satisfy a commutation relation of the form 

[A, B] = a(A)B + b(A) (26) 

for some a(s), b(s) G C[s], then in ( f!^| j it suffices to assume only the n — case, i.e., 

2#o = F(A)^j (27) 

for some F(s) G C[s]. 
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Proof: Since Vo G W, AW C W implies W' := span c {V>o, Aip , A 2 ip , . . .} C W. Since Vo = 
1 + 0(z _1 ) and A raises the order of a function in z by 1, the map W' — > -H+ is bijective, and 
W' E Gr°. In particular, both W and W' are of relative dimension 0, so that W — W. Conversely, 
W = W' clearly implies ipo £ W and AW C W. Assume these equivalent conditions. Then BW C W 
if and only if BW' C W' if and only if the differential equations of the form ( p5| ) are satisfied. Finally, 
when A and -B satisfy a commutation relation of the form (^6|) , the n th equation in ( |25| ) implies the 
(ra + l) st one, so that @ suffices. ■ 

The following propositions take a closer look at the [P, Q] = 1 case and [P, Q] = P case, to show 
that essentially those elements in mJoo which give rise to P and Q in the sense of Lemma |2j and their 
polynomials, are the only elements of the stabilizer. 

Proposition 2 Let p E Z, p > 0. Lei yl E &e suc/i iftai [A, z v \ = 1. J/W £ Gr satisfies 
z p W C W and AW C W, i/ien the stabilizer ofW is generated by z p and A, i.e., S w = C[A, z p \. 

Proof: Since [A, z p ] = 1, A is a first order differential operator in z, so that any C 6 Syv can be 
written as C — 2-oa<i«oo o<j<N a ij z% A^ for some TV > 0. It suffices to prove that = if i < 
or if i mod p. Suppose A raises the order of a function in z by k: ord z Az £ = ^ + k. Let / be 
the set of pairs (i,j) such that i < or i ^ modp, ay ^ 0, and z + kj is maximum among all 
such ay's. We have |/| < oo, and we only need to prove |/| = 0. Suppose this is not true. Let 

°o ■= J2( l ,j)ei CLi o zlA: '- Notm g 

[A.z'AJ] = [A, {z p y' p ]A 1 = {ijp)z l ~ p A\ 

so that ad(A) n (z l A 3 ) = for n ^> if and only if i > and i = modp, we see that for n ^ 
the leading terms of &d(A) n C are ad(A)™Co, which lowers the order of a function in z, and does not 
annihilate the function for a general n. This cannot happen since &d(A) n CW C W, and since in W 
the order of functions in z are bounded from below. I 

Proposition 3 Let p E Z, p > 0. Let A ~ zd/dz — a(z), where a(z) E z + C[[z _1 ]] ; and ipo = 
1 + 0(z^ 1 ) £ 1 + H-. Let W E Gr be the point of the Grassmannian determined by the conditions 
ipo E W and AW C W. Suppose W also satisfies z p W C W. Let F(s) = cIILi( s_c i) e C H> wh ere 
Cj E C, c E C*, 6e £/ie polynomial of degree p as in ( f§^ ) urai/i _B = z p , i.e., V'o satisfies the equation 

F(A)ipo = z p ^ ■ (28) 
Then if F satisfies the following genericity condition: 

( G) For any n ^ mod p, we have (F) + n := XX c » + n ) ^ (-f 1 ) mod p, i.e., n p ((F) + n) 7^ n p (F), 
where (F) = X)f=i( c «) * s ^ e divisor of F, and ir p : C — > C/pZ is i/ie natural projection, 

then the stabilizer ofW is generated by A, z p and £ := z~ p F(A), i.e., 

S w =C[A,z p ,Z}. (29) 

Remark 2 Condition (G) is equivalent to 

(G ') There does not exist n \ p, < n < p, and H(s) E C[.s] of degree n such that F(s) — 
Ut n oH(s-in); 

and if it is not satisfied, i.e., if F(s) — Yii=o H( s ~ * n ) f or some n \ p and H, then taking such 
(n, H) of the smallest n, we observe from our proof below that C[A, z p , £] C Syy C C[^4, z n , where 
I' = z- n H{A). 
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Remark 3 The right-hand side of (p^j equals ■ k>Q Ca l V c k , where [a, 6, c) is any permutation of 
(A, z p ,£); the order does not matter because 

[A,z p }=pz p , [A,£] = -pF(A) and [z p , £] — F(A) — F(A — p) . (30) 
Remark 4 Condition (G) is satisfied by the F in Theorem^: Since 

F ( s )= ( II («-*)) (*-(p-i+cp)) 

and — 1 < c < ; i/iere is no period less than p in the divisor of F modulo p. 

Proof of Prop. [|. Using the commutation relations (^), the definition of W, and Eq. (pq), we 
observe easily that Syy D C[^4,z p ,£]. We prove the converse inclusion in two steps. Only Step 2 
needs Condition (G). 

Step 1. We observe that Sw is spanned by the z-homogeneous elements in Syy, i.e., the elements of 
S w of the form z n f(A), where n eZ and f(s) G C[s]. 

Indeed, let S' C Sw be the subspace of Sw spanned by the z-homogeneous elements, and suppose 
that S" := Sw \ S' ^ 0. Let A*" be a nonnegative integer such that 

gm ■-{ces"\ovd a/az c<N} 

is nonempty. Let C G S"'( Ar ) be such that, writing 

C = 5>V«(4), (31) 

no(C) := max{n | /„ ^ 0} is the smallest in S"( N \ Such a C exists because 
Claim 1 {n (C) \ C G 5"' (7V) } is bounded feeZow. 

Proof: Indeed it is bounded from below by — 2N + 1: since C G S"'W is an ordinary differential 
operator of order < iV, and since -^V'Oj • ■ ■ > -<4-'' v_1 V'o are linearly independent, we have CA l ij) ^ 
for some z, < i < iV. Since A l ij} = (-l)V(l + 0(z -1 )), since CW C W, and since W is a 
span of A J ipo for j > 0, we observe that C does not decrease the order of A 1 ^ in z by more than 
N — 1. This implies, using the notation of (|3l|), that n + deg/„ > — (N — 1) for some n. Hence 
no(C)>n>-deg/ n -(JV-l)>-(2iV-l). ■ 

Now let 

C := [A, C] - n (C)C = J2( n ~ MC))z n f n {A) . (32) 

Clearly C G S'w- We have ordg/a z C" < ord,g/g 2 C < N, and n (C") < n (C) — 1. Hence by the 
minimality of n (C), we must have C S /,( - N \ so that C G 5'. Thus each term (n — n (C)) z™/n (A) 
in ([32]) belongs to 5", and only finitely many f n are non-zero. As a finite linear combination of such, 
we have C" := C - z no{ - c ^ f no{c) {A) G S", so that z"°( c ')/ no(c) (^) = C - C" must also belong to 
Sw, and hence to S', since it is z-homogeneous. This implies C = C" + (C — C") G S", which is a 
contradiction. 

5iep £ Let /(s) ^ be any constant coefficient polynomial, and let n be an integer. We prove that 

z n f(A) G Sw implies p | n , 

and that, when n < 0, z n f(A) G 5w must have the form £ fc /i(A) for fc := —n/p > and some 
/i(a) G C[s]. 

Suppose z n f(A) £ Sw- We assume n/0 without loss of generality. Since z n f(A)ipo G W, by 
Lemma ^ there exists another polynomial <?(s) G C[s], such that 

z n f(A)^ = g(A)^ . (33) 
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First assume n > 0. Let £ > be the least common multiple of p and n. Noting 



z 2 *Vo = z p F(A)ip = F(A - p)z P i> = F(A - p)F{A)4> 

etc., we have 

fl/p-i \ 

[] F(A - ip) 4,0 = z% (34) 
i=0 J 

from (2SS. and 

II G(A-jn) 4x> = z t 4> (35) 
j'=o / 

from (p3]), where G(s) = g(s)/ f(s — n) is a rational function in s, and G(A — jn) in ( |35| ) is understood 
as an element of the field of fractions of C[A]; this makes sense because, since {A n ipo}n=o,i,... 1S linearly 
independent, the representation 

CM 9 m ^ /(A)Vo e W 

is faithful. 

Comparing the left-hand sides of j34| ) and (^), we thus have the equality 

i/p-l l/n—l 

f[ F(s-ip)= G{s-jn) (36) 

i=0 j=0 

of rational functions in s. Since the left-hand side of it is a polynomial of s, so is the right-hand side. 
Let D be the divisor of this polynomial, and let ire be the natural map C — * C/£Z. From the left- 
(resp. right-)hand side of (^) the image ire(D) of divisor D on the cylinder C/£Z is invariant under 
the translation by p (resp. n). But the genericity condition (G) implies that if tti(D) is invariant 
under the translation by k £ Z, then p \ k. Hence p \ n. 

Note here that, since £ is the least common multiple of p and n, this implies £ = n, so that the 
right-hand side of ([56]) is G(s) itself. Hence 

n/p—l 

g(s)/f(s-n) = G(s)= f[ F(s - ip) . 



In particular, g{s)/ f(s — n) is a polynomial. 

In the case where n < 0, after rewriting (133) as 



i=0 



z-" ff (A)Vo = f(A)4> , 

we switch the roles of / and g, and n and — n, to proceed exactly the same way to prove p \ n and 

-n/p-i 

f(a)/g(s + n)= ]J F(s-ip). 

i=0 

Thus we have 

/-n/p-l \ 

= W n F(A-ip)\g{A + n) 

= (z- p F(A)y n/p g(A + n) 
= Z k g(A + n)=:t k h(A), 

proving the last assertion of Step 2, and hence completing the proof of Prop. |3|. 
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2.2 Symmetric functions and matrix integrals 

In this subsection, we prove a number of lemmas regarding symmetric functions. 

Lemma 4 Let s and N be positive integers. Let F(x^\ . . . , x^) be a function which is symmetric 
in each x^ :— (x^, . . . , x^ ) <E C^, r = 1, . . . , s; let fx, . . . , f s be functions of two variables, and let 
B(x^) be a skew- symmetric function of x^ s \ If Ci, . . . , C s denote s fixed contours in C, then the 
integral 

s N 

(r) 



:=/•••/ nn^ 

J i(C 1 )»x...x(C s )\ =l!=1 

S 

■F(x {1 \ . . .,x {s) )B(x^) Y[ det (/,(.. . . ,< 



(r-i) Jr)\ 

J ' J l<i,j<N 



where x^ € C N comes in as the first argument of fx, is skew- symmetric in x^°\ and 



<5>(x^) = (Nl) s f ■■■ f Tidx^ 
J J(C 1 )«x--x(C s )« ^7 

s N 

■f(x^, . . „,w)b(xW) n n /rK (r " i u (r) ^ 

r— 1 i—1 



Proof: For any (good) functions A = A{x^\ . . . , x^) and h — h(x^ l \ . . . , x^), let 



(Ah);= I ■■■ I T\dxP>A(xV\...,xM)h(xW,...,xW 

'(Ci)«x---x(C s )« *| 



For any o> € S N , let x$ := (x^, . . . .arj.^), and (x^, . . . , x^) := h(x$ ,...,x£). 

Clearly (Ah) = /A^" 1 ''"' (7s 'hS' ri '"''' 7s ^'). If, moreover, A is symmetric in each of x^\ . . . jX^ 1 ), and 
skew-symmetric in x^ s \ i.e., A^" 1 '" — (— ly^^A, then we have 



(Ah) = (^•-'"•W* 1 --"')) = (_ 1 ) E K)^ /l ( CT i,..,^)^ > V(Tj , e e 



JV ■ 



Applying this to h(x^\ . . . , x^) :— ]J r ]J i f r (x\ r ^,2:^), and summing it up over (a±, . . . , a s ) £ 
(&n) s , we obtain 



(NiriAinifri^,^) 

r i ' 

A ^ (-l^nn/^o^)- witha =id, 



Setting here A — F(x^\ . . . , x^)B(x^) proves the identity in Lemma ||. Finally, ^>(x^) is skew- 
symmetric in x^ since det (fi(x\°\x^)^ is. I 
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Lemma 5 (See ^ Lemma 4.2], @ Eq. (2.21)], /§§ Theorem 8.18].) Let 

W = spa,n c {ip (z),ip 1 (z),il)2(z), ...} € Gr 

with functions 

ipk(z)= X] a 3,k zJ > fc = 0,1,2,..., 

— oo<j</c 

smc/i that atk = 1 /or k ^> 0, i.e., ord z ipk{z) < k, and ipk(z) — z h (l + 0(z^ 1 )) for > 0. Let N > 
be any integer such that this condition holds for k > N. Let z\, 
near oo. Then the t -function r(t) at 



z~ n 



zn be formal scalar variables 

(37) 



is given by 



det^-i^)) <JV 

T(t) = — =-!^=— 

det(z? )i<%,j<N 



(38) 



Proof: Our proof is based on Kontsevich's idea in [flSfj ; see [[i~7| Sect. 2.3] for a proof using free fermions. 
To keep the notation simple, let us denote by (1 — z)^ 1 and (— z + 1) _1 the geometric series z n 
and — X)-L z ™' respectively. Let 5(z) := (1 — z)^ 1 — (— z + 1) _1 = X^oo z ™' which plays the role of 
delta function, in the sense that 

5{z/y)f(z) = S(z/y)f(y) , (39) 

as is obvious by taking f(z) = z m (see Let a := Y\f =1 (— Zj) = (—1) n Z\...Zn. Let Oi := 

— z i/ z j)j i = understood as rational functions of Zj's, so that we have the 

following identity of formal power series in z: 



N 



N 



From (|37|) we have 

3 := exp 



n=l 



JV 



i=l 
iV 



i=l 
N 



N 



i=l 



N 



i=l 
N 



5>*(*/*i)+n 





0- 






z. 


(- 


— -+ 











+ 1 


-1 

1 









so that by using (p9|), we have 



i=i Vi=i ^ Zj 



i=l 



= ^ Vitiz/ziWjizi) + z- N (a + 0(z- x ))^{z) . 

i=l 

Denoting by B the matrix of the composite map in (|Tj) with respect to the bases {V^'lj^o an( i 
{z k }kLo' we nave tnus B = B° + B l , where 



/ 1 



5° = 



6 i 

..-2 



N 



'JV 



/ ipo(zi) ipx(zi) ■■■ 
\ V'o(zjv) ^i(zn) ■ ■ ■ 
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JV 



B 1 








a 

a 



a-2,o 


a-2,i ■■■ 


a-i,o 


a-i,i ■■■ 




Ooi 





an 









Sn is the diagonal matrix diag(cri, . . . ,(?n), and akj, — oo < fc < oo, < j < oo, are the Laurent 
coefBcients of ipj = ^2 k akjZ k . 

Let us apply some column operations on B. Adding an appropriate linear combination of first N 
columns to the (N + i) th column (i > 0), we can eliminate the column t (^pN+i(zi), ■ ■ ■ , ^Pn+i{zn)) , 
i > 0, from B°. Since N is large enough so that ajj = 1 for j > N, B 1 has the form 



a 



ooxjV 



V 



a * \ 

a 







/ 



so that the "*" part can be eliminated by further column operations on columns N + 1, iV + 2, 
which do not alter the £>°-part. Here O mxn is the m x n zero matrix. The matrix B can thus be 
reduced to B' = B'° + B' 1 , where 



B 



/o 



V : 




Nxoo 



Let n, n > TV, be an integer. Note that the column operations needed to bring B into B' only adds 
linear combinations of lower numbered columns to higher ones. Hence, denoting by B n , B' n , B'® 
and B'n the matrices of the first n rows and columns in B, B', B'° and B' 1 , respectively, we have 
detB„ = det B' n = det{B'£ + B£), with 



/ 1 



B" 



-n+l 



-n+l 




iPn(zi) 



o 



JVx(ra-JV) 



/ 



and 



BL 



0( n -N)xN 


aln-N 


OnxN 


Onx(u~N) 



Since the last n — N columns of B'® arc 0, we have 



BL 



crI n -N 



O 



Nx(n-N) 



where Z consists of the last N rows and the first N columns of £?'°: 



( H 



■n+N 



-n+N 



z = 



N 



n+l 



-n+l 



JV 




4>n{zi) \ 
-0Jv(^jv) / 
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Hence we have, using a = {—l) N Zi . . . zn, 

dctB n =detB' n = (-1)^"-^) det Z det(aI n - N ) 
= { Zl ...z N ) n - N detZ 



( Zl ... z N f- ' N det Z' , 



where 



(4 



N-l 



z' = 



y N-l 
'N 



'-N 



s 



N 



I ^o(zi) ■■■ iPn{zi) 
\ iPo(zn) ■ ■ ■ iPn(z n ) 



Noticing 
and 



det S 



N 



\i •■■ i / 

det^)^,^ = (-l) Ar(Ar - 1)/2 det(zf 1 ) 1 < lJ < JV 



JV-l 



N 

i 



[zi... z N 



,N-1 



). Since n > N is arbitrary, this 



we observe that det B n coincides with the right-hand side of (|3_ 

completes the proof of Lemma |5[ I 

Lemma 6 Let Z := diagOi, . . .,z N ). Let A := ((p - 1){N - 1), (jp - l)(N - 2), . . . ,p - l). For a 
polynomial f(y,z), let us denote by (/(y, z)) 2 the terms in f(y,z) which are quadratic in y. Then 
we hav^ 



A(zP) 



Fx f -trZ, tr Z 2 . tr Z 3 , 
'2 '3 



n 



dF exp tr — 



"Hat 



p + i 



where c is a non-zero constant which depends only on N and p. 
Proof: The Schur function associated with the partition A is given by (see 

. A A+5 (y) 



(TV JV ^ Af 

1 1 1 



A 5 (y) 



where 5 = (iV — 1>JV — 2>>-<>l>0) and A M (y) = det(yf J )i<i,j<iv- Therefore we have, with 
A + 5 = (p{N - 1) > p(N - 2) > • • • > p > 0), 

A( z p) A x+5 (z) 



F x 



/ N 1 N l N \ 

(-E*'-2E*i'-3E*?""J' 



A(z) A«(z) 

establishing the first equality of Lemma |6j In order to establish the second one, note 



tr 



HY + Z)p+ 1 
V P + l 



1 



= -trfy 2 ^" 1 
2 1 



yzyz p - 2 + --- + yz p - 1 y) 



i<\ is the Schur function for the partition A. 
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Hence, performing a Gaussian integration, we find 

^ex P trf- (r + Z f +1 ) = fdYeJ-Wy^ Zf ~ Z ' 



l'-'r I ./ " "* V 2^r!' Xijiji Zi-Zj 



(2-) Ar2/2 n 



1/2 



z p - z p 

\<i,j<N 1 3 , 



(2^) jV2/2 TT g - gj A 
n N/2 11 7 P _ 7 P 1 1 

P l<i<j<N 1 J 1 

(2<^A(z)n 

pAT/2 A ( z p) 11 < 



2„- 



establishing Lemma 

Remark 5 In general we have 



dY e -HV(Y+z)) 2 =(27r )^ 2 /2 A (*) 



A F 



The following lemma is due to Harish Chandra, Bessis-Itzykson-Zuber and Duistermaat-Heckman 
among others: 

Lemma 7 Given N x N -diagonal matrices X and Y , we have 

tr xuyu^tj ( „^im=2l det(e^^)i< M < A 

A proof can be found in . 

3 Matrix Fourier Transforms 

In this section we explain how generalized Kontsevich integrals (see (l9[ [|, |24j ) are closely related to 
the theory of Fourier transforms. Indeed, if V(x) grows sufficiently at infinity, any Fourier transform 

/>oo 

a(y) = / e- vi - x)+xy dx (40) 



leads to a linear space of functions W invariant under two operators A and V'{z) satisfying [A, V'(z)] = 
1. 



V'(^)a{ V ) = ,,„(,,). (41! 



(i) The point is that a(y) satisfies the differential equation 

' d_ 
dy 

as seen from 

/CO o poo 
°-v( x)+X y dx = / (_v"(a;) + y)e- v ^ +xy dx 
-co J -co 

Thus setting y = V'(z) in (|l|) and A := V'^z^d/dz = d/dy\ y=v , (zy the function a(V(z)) 
satisfies the differential equation 

V'(A )a(V'(z)) = V\z)a(y'{z)) . (42) 
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(ii) The method of stationary phase applied to integrals (40) and their derivatives leads to the 
following estimate, upon Taylor expanding V(x) around x = z, 



f-i \ n 
O \ 

By 



x n e -V(x)+xV'(z) dx 

x n e -(V(z) + (x-z)V'(z) + (l/2)(x~z) 2 V"(z)+0(x- Z ) 3 )+xV'(z) dx 



= e -V(z)+zV'(z) / x n e -(l/2)(x-z) 2 V"(z)(l + (V"'/V")0(x-z)) dx 



-V(z)+zV'(z) _ 



'V 



(L(^ +z f c "'' 2d « +o{i/z) 



p{z)- l z"(l + 0(l/z)). 



with 



Therefore defining 



and 



P(z) 



^V(z)~zV' (z) /T7JJ 



p(z)- 



V=V'(z) 



ip n {z) := A n ijj (z) := p(z)-^a{y) 



n = 0,l,. 



y=V'{z) 



the differential equation ( |42|) implies 

V'(A)Mz) = V'(z)Mz)- 
This, combined with (|43|), proves that the linear span 

W := span c {^(z) = z k (l + 0(l/z)) | k = 0, 1, 2, . . .} 
is invariant under the operators A and V'(z), i.e., 

AW <ZW and V'(z)WcW, with [A,V'{z)) = l. 
(iii) By Lemma |5|, the r-function corresponding to W, at time t as in ([37]), is given by 



r(t) = 



det(^ )i<i,j<jv 



A(z) 



det p(z; 



j'-i 



e _y( a: )+xj/ da . 



j/=v"( Z «) y 
'(*«) 



l<iJ<AT 



AHA(z) 



l<a,0<N ' 



using Lemma ^ with s = 1 and the skew-symmetry of A (a 



jV!A(z)/A(F'(z)) V 



detfe 3 ^'^) 



l<q,/3<jV 



A(ar)A(y(z)) 



(43) 
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A(z)/A(V'(z)) V 1 ; 7u ( .v) 

using Lemma 0, with X = diag(x), 

- c ' e tr(v(z)-zy'(z))Jw« I I usine Lemma S 

ce f dXe _tr(y(x+z)) 2 ' usmg i^emma |C|, 

= C " f ~ d y e -tr(v(y+z)) a ' upon setting X = r + 
for some constants c, c' and c" depending on N . 

4 Generalized Hankel Functions, Differential Equations and 
Laplace Transforms 

This section deals with the properties of Hankel functions and their generalizations. 
Lemma 8 The family of integrals 

c+l roc -c( uz \k e -{u-l)z 

admits, for large z > 0, an asymptotic expansion in C((z -1 )) o/ i/ie /orm 

V> fc (z) = z fc (l + 0(l/z)), (45) 
with ipoi z ) satisfying the differential equation 

fV ~ l ' d \ ^ / 8 



. dz J \ dz 

si=0 v 7 i=0 



eZz " c ( II ( Z 7T, - 1 ) - cp II ( - ) zC e~ z Mz) = {-z) p Mz) , 

or equivalently 

e z z~ c ^z* (J-J - cpzP- 1 (J-J j z c e- z Mz) = (-z) P Mz) . 
Moreover ipk{ z ) admits the following representation in terms of a double integra^ 

„c+l r roc 

= z (p-i)(c+i) e z / dajx^-^ / dy f k (y) e - xyP , (46) 

2 ™ Jo Jo 

where, in the first integral, 7 := 7+ + 7~ C C denotes the contour consisting of two half-lines 
7 ± = R + £ ±:L , C := e* 1 /' 9 ', through the origin making an angle ±7r/p with the positive real axis, with 
the orientation given as to go from C^ 1 • 00 to to £ • 00 (see Fig. [J (a)), and where in the second 
integral, 

My) = (c k+1 e~<y - r*- x e- rl V = E 

3=0 J ' 

where a n = C™ — C~™ = 2isin(n7r/p). 
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(a) 



(b) 



Figure 1: Contours of integration: (a) contour 7 for (f46|); (b) closed contour for jig) 



00 e - v 



Proof: Setting v — (u — l)z, and using 

r(-c) = y —dv for C <o, 

we first observe that for each n > 0, 

, , , p c+1 z fc f°° (l + v/z) k e- v 
Vk(z) = r/ ~s / ; 77+T dv 



r(_c) y ° ^ C+ V +1 (1 + 1 (SK*/*)''" 1 ))' 

= z fe (l + fc^z- 1 + • • • + 6 fc ,„^- n + 0(l/z n+1 )) 

as z — > 00, where the bk,i '■= (T(—c + i)/T(—c))bk,i — (jTj=o(~ c + j)j &k,» are obtained from the 
coefficients bk,i of the expansion^ 

" +S,A r = l + f>,^, (47) 



i + i(El 2 (I)^)) c+ 



confirming the asymptotic expansion 
Moreover, setting 



we have for c < and Re z > 

= -z p ~' 



(up - 1) 

, 0u V ~'(uP-l) c+1 



«=i 



5 If p = 2, so that 7 becomes the imaginary axis, these integrals should be interpreted by replacing £ by £ e = 
e (7ri/2)- £j and 7 by + R+f- 1 , and then taking the limit as e I 0. 

6 Noting that the radius of convergence of this power series is | £ — 1 1 , one can get a precise growth estimate of the 
coefficients of ipk( z ) which implies that, in particular, as always with the string equation, W does not belong to the 
L 2 -Grassmannian of Segal- Wilson 
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-c^ — uz 



(-If I ({-zuf cpi-zuf- 1 - {-zY) { * p \ )c+1 du 



\i=0 

using in the last line the operator identity 



n(*|-')-"(s' 

i=0 v / \ 

thus showing that V'o( 2; ) satisfies the differential equation (||) or (j^'). 

Consider a bounded domain DcC, whose boundary consists of the lines 7^, making an angle 
±n/p with the positive real axis, two circle segments C fi , about the origin, of large enough radius R 
and a small circle about 1 of radius e connected to , as in Fig. 0(b). The function e - uz /( U P- l) c+1 
is univalued in D and all its singularities lie outside D. By Cauchy's theorem we have 

r r R r \ (uz) k e - ( - u - 1 >- 



(uP - l) c 



du = . (48) 



h R J~t + R JC R JR JC E Jl+e JC~ 

Observe that, for z > and p > 2, we have zcos6> > zcos(ir/p) > for < 8 < ir/p, implying 
[ ± ^ e ~ ( "~+" du = (i^-(c+i)P+i e -fl-os(Vp)) ^ o 

as R | 00. Since c < 0, we also have 

11 du = 0(e~ c ) ^0 



as e I 0. So, taking limits as e J. and i? | 00 leads to 



(up-i) c+1 V7oo A-io / (« p -i) c+1 

= (^^(c+i)) / e_> 

(uz) k ; P -("- 1 ) z 
= 2ze-" c sin7rc / 1 J - - ^ . 



(up - l) c +! 

Note that, since u p — 1 < along 7, we have the following T-function representation 



and thus 



P 



dxx c e x ^ P - 1) 



r(-c) y x (uf - i) c +! 

p c+l e 7ric _ c /■ ( uz )fc e -(u-l)2 



-Z 



du 



2ismircT(-c) 7 7 (nP - l) c+1 



p c+i z -c 



2isin7rcr(-c)r(c+ 1) 



^~z- c I du {uz) k e- {u ~ 1)z / </.,■ ,-' , ' ; 

____s(p-i)(c+i) I dw I dxe z - w w k x c e xw "e~ xzP 
2m 
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upon setting w — uz. Here we used the T-function duplication, T(— c)T(c + 1) = — 7r/ sin7rc, —1 < 
c < 0. Working out the integral over 7, interchanging the integrations and using we find 

n c+l roc 

^Jz) = t z (p-V(c+l) e * / dxx c e- xz " ■ 



p0 poc 

k ^ dye- rly y k e- xy " +( k+1 dy y h e- xyi 

Joe JO 



ith 



P_ z (p-V(c+l) e z J dxx c e~ xz " J -/.// /,.(/;)( " 



f k (y) = (e +1 e- iv -C k - 1 e- rly )y k 



as announced in ( f46| ) , thus ending the proof of Lemma ^. 

Lemma 9 The linear space spanned by the generalized Hdnkel functions, 

f P c+1 f°° Z~ c (uz) k e~^ u ~ 1 ^ z 

W = span c |Vfc(*) = Y(ZcjJ i {uP _ 1)c +i du & = 0,1,2,. 

is invariant under 

d d 
z p and A c := z c e z z— o e z z c — z— z + c 

fso i/iai z p ] = z p ), with tpQ satisfying the differential equation (Qj. 

Proof: The space W is invariant under A c , because 

. . . , P c+1 _ c z 8 c _ z f°° z- c (uz) k e-( u -V z , 

- ^ z-e'zl r du 



Moreover, the operator 



r(-c) dz J 1 (up - l) c +! 

ki> k {z) - ip k+1 



p— 1 p-2 

l[(A c -i)-cp]l(A c -i) 

i=0 i=0 



has the form X)o a i^c; with ct p — 1. From Lemma |[ the solution to the differential equation 

/p-i p-2 \ 

[l[(A c - i) - c P ]l(A c - i) i> (z) = {-zfMz) 

\i=0 i=0 / 

is given by the function in ( f44| ) or ( |46| ) for = 0. An asymptotic expansion of the form 

^ O (z) = l + O(z- 1 ) 

follows from (E5j). 
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5 Proof of the Main Statements 

5.1 Proof of Theorems |3| and [I] and Remark [I] 

In Lemma |^, we have constructed a space W and an operator A = A c such that 

AW C W and z p W C W , 

with the lowest order element -00 S VV satisfying Eq. (@) . Proposition || and Remark || imply that 
the stabilizer of W is CL4, z p , z _p P(A)], proving Theorem |, Part (i). 

Let \1/ and r be the wave function and the r-function, respectively, associated with the KP time 
evolution W* = e^ f - 2 'W of W. We now define the operators P and Q in the x-variable, via the 
operators A and z p in the z-variable, by means of 

z p -$(t,z) =PV(t,z) and (l/p)AV(t,z) = QV{t,z). 

According to Lemma|, P and Q are differential operators. They satisfy [P, Q] = P since [(l/p)A, z p ] = 
z p . Note that P and Q can also be written: 

P = L p = SDPS- 1 

and 

Q = - (ML - L + c) = -S | V - £> + c ) ST 1 , 

P P \ i / 

where 

r(t-[^]) 

in terms of the r- function above, and L and M are as in ( |l2| ) and (|l6|), proving Theorem Part (ii). 

Since (M — 1)L = pQ — c is a differential operator, we also have, using the notation a^- as in the 
statement of Theorem |3|, 

((M - l)L) m L np = a m,i( M - iy i L" +np 



0<i<m 
0<j<i 



c 

is a differential operator. Thus 



J2 a m,i[ .)(-iy- j (MiL i+n v)^ = 0, 



0<i<m 
0<j<i 

implying ([To]), upon using ([l9|), completing the proof of Theorem ||. 
To prove Remark [l], we evaluate 

at t = to find — ((pT / dt\) / t) | t = a. Remember, on the one hand, 

4<(0,0,z) = 7/;o(z) = (1 + 601^ + ■■■)- 

and on the other hand 



*(0,0,s) 



r(ii + x — z 1 , . . .) 



r(ti +x,. ..) 
-1 <5 T -1 



a;=0,t=0 
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Therefore a — boi = (— c)&o,i = c(l + c)(p — l)/2 as stated in Remark [I], as implied by ([47|). 
To prove Theorem [j], note that at t = 0, 

Q| t=0 = (l/p)S((a;-l)(5/fe) + c)5- 1 

= (l/p)(x — l){d/dx) + c + (negative order terms) . 

Since Q must be a differential operator, the negative order terms vanish, and Q\t—o — (l/p)(x — 
l)(d/dx) + c. Thus, from the second equation in (^), we have 

0= (Q\ t=0 - A c )<i>(x,0, z) = ((x - l)(d/dx) ~ z{d/dz - l))^{x,0, z) . (49) 

Since this is a first order equation and the line x — is noncharacteristic, ^(x,0, z) is determined 
by (|49| ) together with the initial condition $(0, 0, z) = ipo(z). It is easy to check that the right-hand 
side of (||) satisfies these conditions. Finally, (Q) follows from (||): Writing (@) as 



substituting (1 — x)z for z, using the scaling invariance of zd/dz, and dividing both sides by z p , we 
get 

(50) 



^f(z^ + c) (e^'MO- ~ x)z)) =(x- lye^-V'MQ- ~ x » ■ 



Multiplying both sides of this formula by e z , and using the identity e z (zd/dz+c) = (z d/dz—z+c)oe z , 
we get the second formula in (0). Next, switching the roles of z and 1 — x in (^0j), we get 

^_L-f((x - 1) A + ^ ( e (*-D^ o((1 - *)*)) - feb-V'MQ- ~ x)z) ■ 

Multiplying both sides of this formula by e z , and using the fact that e z commutes with (x— l)d/dx+c, 
we get the first formula in (|?|) , completing the proof of Theorem [j]. 

5.2 Proof of Theorem [2| 

Setting t n = —— n = 1, 2, . . . , and using Lemma |!|, and Lemma ^ with s = 2, we have 

det (^-i(zi)) 1<k N 



r(t) 



A(z) 

N / /-co />oo \ 



A(z) 

a N S 2 (t) 
A(z) 

a"S 2 (t) 



/ dy(r\xA 



(iV!) 2 A(z) 

■det(/ fc _ 1 ( W )) M det (e" x ^ ? ). .det (e"^) 
a Ar S' 2 (i)A(^ 



4*4* (n«) 



(iV!) 2 A(z) 

det (e- 3 ^) det (e^ 1 ^) 



A(x)^A(y)^ 



■So(y) 



A(x)A(zp) A(x)A(yf) ' 
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where a = p c+1 /2iri, 

N 



1 

and 

A /-.r>\ HptV f, , i 11 V, 

>l<i,k<N 



A(yf)det(/ fc _ 1 (y l )) ] 

bo{yi,V2, ■■■,vn) = 



My) My) 

So we have, for some constants C, C and C" depending on N, p and c, 

r(t) = cr frWy) f d:E A(x) 2 / dyAM'Soiy)- 
A(z) Jr™ im™ 

U(JV) Ju(N) 



using Lemma [7] 



c S 2 (t)A(z^ " 



U(N) JU(N) 



setting Uy — VyUx for fixed Ux in the last 
integral and noting that dUx AUy = dUx dVy 

= c S2{t } A [ zP) [ dxA(xf(f\xt) [ dU x e-^?*v x 



A(z) 



U(JV) 



dyA 2 (y)S (y) / dU Y e~ uu - xUxU ^ ^ 



c ,S 2 (t)A(zP 



C"5i(t) 



/ dX det(X c )e~ tl ZPX [ dYS (y) 
A ( z ) •/«+ 

/ w+ dXdet(X c )e- trZPX / w+ drS (y)e- trXyP 



e -trXY» 



J n dX exp tr 



•Hw * V p +1 



where we used Lemma ^ in the last equality, and the definition of Si (t) in Theorem ^|. A similar 
calculation, outlined below, implies the second formula for r, upon using the first representation of 
rl> k (z) in ©: 

detU*- 1 *^,*))^, 
T W = ~A(z) ~ ~ ' " = ~" ^ i=1 Zl ' 



A( 

N 

« e E „TT z ( P -i)(c+i). 
A (2 



^-det fae^4 p - 1)(c+1) / dw [°° dxe^w^ 1 

(z) V A Jo 

n 

/• /• />0O />OC JV 

■■■ dw ■■■ dx e~ E ^ TT x ■ AH e - *« 184 TT e : 
J 7 J 7 Jo Jo i=1 j=1 

QjV e I>»TTz (p ~ 1)(c+1) — /" / dxe-^ Wi T\x?A( 



k.i 



■det(e- ft )Ku<« det (^ 



l<i,j<N — V Jl<i,j<N 
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f^j m(dW)J n+ dXdetX c (A{ w P)/A(w))e t ^ z - w h t " x ^ WP - z ^ 



Jj. dX exp tr 



(X + Z)P+ 1 

p+i /2 



ending the proof of Theorem |^. 
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